Supersymmetry and Goldstino-like Mode in Bose-Fermi Mixtures 
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Supersymmetry is assumed to be a basic symmetry of the world in many high energy theories, but 
none of the super partners of any known elementary particle has been observed yet. We argue that 
supersymmetry can also be realized and studied in ultracold atomic systems with a mixture of bosons 
and fermions, with properly tuned interactions and single particle dispersion. We further show 
that in such non-releativistic systems supersymmetry is either spontaneously broken, or explicitly 
broken by a chemical potential difference between the bosons and fermions. In both cases the system 
supports a sharp fermionic collective mode similar to the Goldstino mode in high-energy physics, 
due to supersymmetry. We also discuss possible ways to detect this mode experimentally. 



Supersymmetry, which is a symmetry that relates 
bosons and fermions, became a topic of strong interest 
in elementary particle physics after Wess and Zumino 
constructed the first "realistic" model f^. It may play a 
fundamental role because the supersymmetry algebra is 
the only graded Lie algebra of the symmetries of the S- 
matrix consistent with relativistic quantum field theory 
[2] . The supersymmetric string theory is a unique theory 
expected to give a unified description of all interactions 
in nature J^. However, none of the super partners (which 
have identical properties except for opposite statistics) of 
any known elementary particles has been found in exper- 
iments thus far. Therefore, it is extremely important to 
study the breaking of the supersymmetry. 

Recent experimental progress in mixtures of ultracold 
Bose and Fermi atoms '4| provides an opportunity to re- 
alize and study supersymmetry in such atomic systems. 
Theoretically, several different models that exhibit su- 
persymmetry have been studied. An ultracold super- 
string model was constructed The physical behavior 
of an exactly soluable model of one-dimensional Bose- 
Fermi mixture was investigated 0] . A general formalism 
to study such supersymmetric systems based on coher- 
ent state path integral was set up in Ref. 0. In a recent 
work Q, we studied a supersymmetric Hubbard model 
and focused on the Mott insulator phase for bosons. 

In the present paper we study some general properties 
of supersymmetric Bose-Fermi mixtures, in which bosons 
and fermions are supersymmetric partners of each other. 
We show that in the presence of time-reversal symmetry, 
supersymmetry is always broken, either spontaneously 
or by a chemical potential difference between bosons and 
fermions. For a truly supersymmetric grand canonical 
Hamiltonian, we find the system contains bosons only 
even though the grand canonical Hamiltonian is invari- 
ant under a supersymmetry transformation that turns a 
boson into a fermion; in this case supersymmetry is spon- 
taneously broken. To support a finite density of fermions, 
one needs a chemical potential difference between the 
bosons and fermions: A/i = /ii? — /iB>0, which breaks 
supersymmetry of grand canonical Hamiltonian explic- 



itly but keeps the canonical Hamiltonian supersymmet- 
ric. We find in both cases the system supports a sharp 
fermionic collective excitation similar to the Goldstino 
mode in supersymmetric high energy theories, which is 
gapless in the former case while has a gap equal to A/i 
for the latter. For simplicity we refer to this Goldstino- 
like mode as "Goldstino" in the following. We will also 
discuss its possible experimental detection. 

While our results are general, in the following wc illus- 
trate them by considering a simple lattice model (in the 
grand canonical ensemble) with mixture between a single 
species of bosons and a single species of fermions: 



H = f + V; 



(1) 
(2) 
(3) 

(4) 



Here and fi are the boson and fermion operators on 
site i, and nf and n{ are the corresponding number oper- 
ators; Nb = '^i' = J2i presence of 
time-reversal symmetry, all the hopping matrix elements 
(t's) are real; we further assume they are all non-negative 
so the hoppings are not frustrated as is usually the case. 
We now introduce generators of supersymmetry Q: 



(5) 



It is easy to verify that they satisfy the following rela- 
tions: Q2 ^ (Qt)2 ^ 0; = N ^ Nb + Nf; and 
[Q, N] = [Q^ , N] ~ 0. From these relation it is clear that 
Q is a fermionic operator. Physically Q turns a fermion 
into a boson, and does the opposite. While in (l5|) Q is 
defined in 2nd quantized notation, we also need to know 
its operation on a 1st quantized wave function 



■ ,XAr_g;yi, 



(6) 



where x and y are the coordinates of the bosons and 
fermions respectively, so that ^ is symmetric under the 
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exchange of x's while antisymmetric under the exchange 
of y's. We find (up to a normalization constant) 

(3V(xi, • • • jYa'f) = 5'xi,...,x„g,yiV'(xi,---,ywF)>(7) 

(3'l'^(xi, • • • ,yAr^) = Axi,yi,...,y„^-0(xi,---,y7VF),(8) 

where 5'xi,- -,xjv3.yi is the symmetrization operator for co- 
ordinates Xi,- • - jXjv^jyi, and ^xi,yi,- -,yNp is the anti- 
symmetrization operator for coordinates xi, yi, • • • , yjvj? . 
Physically, S turns a fermion into a boson by symmetriz- 
ing one fermion coordinate with respect to all boson co- 
ordinates, and A turns a boson into a fermion by anti- 
symmetrizing one boson coordinate with respect to all 
fermion coordinates. 

If the bosons and fcrmions have the same disper- 
sion and interaction strength, namely tfj = tfj, and 
it is easy to show that the canoni- 
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cal ensemble Hamiltonian H is supersymmetric, i.e. 
[Q,H] = [QlH]^0. 



(9) 



For completeness we also present a generic example of 
supersymmetric H in the continuum, in 1st quantization: 

i<NB i<NF 

i<j<NB i<j<NF 

+ J2 U{^.-y,), (10) 

i<NB,j<NF 

where V is single particle potential and U is two-body 
interaction. Of course the 2nd quantized, lattice Hamil- 
tonian (2) also has a corresponding 1st quantized form: 



Hist — Tist + Vist , 



(11) 



where the form of Ti^f and Vist can be deduced from Eqs. 
(3) and (4). If ?/'(xi, • • • , y^v^ ) is an eigen wave function of 
the Hamiltonian Hist with Nb bosons and Np fermions, 
then (5'0(xi, • • • jYNf) is an eigen wave functionof Hist 
with Nb + 1 bosons and Np — 1 fermions with exactly 
the same eigen energy. Similarly Q^tp{xi, - ■ ■ ,ypfp) is 
an eigen wave function of Hist with Nb — 1 bosons and 
Np + 1 fermions, also with exactly the same eigen energy. 

If we further have the same chemical potential for the 
bosons and fermions: fip — fiB, the grand canonical 
Hamiltonian is also supersymmetric: 



[Q,Hg] = [QIHg]^0. 



(12) 



Given the great tunability of parameters in cold atom 
systems, we expect such conditions can be reached in a 
variety of systems. In the following we discuss conse- 
quences of such supersymmetry when present. 

We start by considering the case where Hq is su- 
persymmetric: [Q,Hq] = 0. In this case we can 



prove that the ground state of Hq contains no or only 
one fermion. Let us assume the ground state con- 
tains more than one fermion, and has the wave function 
'0(xi, • • • ,xjvB;yi, • • • ,yjvF)- Due to the time reversal 
symmetry of Hq, ip can be chosen to be real. Now con- 
struct a trial state: 

'0(xi,---,XAr^;yi,---,yAr^) = |'0(xi,- • •,yjv_p)|. (13) 

Obviously, tp is non-negative and different from ij: for 
Np > 1, as the latter changes sign under the exchange 
of y's. In fact V' is a two-component or "spin"-l/2 bo- 
son wave function as it is also symmetric under the ex- 
change of y's. Since V' and ij} differ in phase only, we 
find the potential/interaction energy does not change: 
(V'lVistIV') = {4'\Vist\4') , because Vist depends only on 
density but not the phase of wave function. The situa- 
tion is different for Tist, which is sensitive to the wave 
function phase. Because tp and ip are the same for certain 
configurations (when ip > Qi) and differ by a — sign for 
others (when -0 < 0), we find the expectation value of 
every term in (Tist) is either the same for {ip) and lip), or 
differ by a — sign. Because ^p is non-negative, and all t's 
in T are non-negative (meaning T only has negative ma- 
trix elements), we find every single term in {'p\Tist\'p) is 
non-positive. We thus have {'ip\Tist\'p) > {ip\Tist\ip) , and 
as a result {tp\Hist\ip) > (iplHistltp) ■ Since in general I'p) 
is not an eigen state of Hist , and the ground state of such 
spin-1 /2 boson system is spin fully polarized 3, 11 1 , from 
variational theorem we find (i/jlHclip) > Eq, where Eq is 
the lowest Hq eigen value for the case with N bosons and 
no fermion (here we also used the fact that bosons and 
fermions have the same chemical potential). This is in 
contradiction with the assumption that \ip) is the ground 
state of Hg- We thus conclude that in the ground state 
Np can only be or l[l^. 

We now show that the ground states have a double 
degeneracy. Assume \tpo) is the ground state with Np = 
0. We can then construct a different state with Np = 1 
and one fewer boson 
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IV-o), 



(14) 



which is also an exact eigen state of Hq with exactly 
the same energy Eq, as guaranteed by supersymmetry 
(1121) . IV'g) can be viewed as a zero momentum, fermionic 
zero mode of the ground state, which is known as the 
Goldstino mode in the high-energy literature. 

As a simple example illustrating the results presented 
above, consider the special case of non-interacting par- 
ticles with all U — 0. For non-interacting bosons we 
always have ^b — (we measure energy from the bot- 
tom of single particle dispersion) at zero temperature, 
regardless of boson number. Supersymmetry of Hq re- 
quires fJ-F = fJ-B =0, as a result we can have either no 
fermion, or a single fermion occupying the k = state 
(with zero energy) in the ground state of Hq- 
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It should be clear from the discussion above that not 
only the ground states, but all eigen states of He (ex- 
cept for vacuum with = 0) come in degenerate pairs 
that have the same total number of particles but differ 
in fermion number by one due to (|12p. The vacuum is 
special in that it is the only state that is supersymmetric, 
as it is annihilated by both Q and . 

We can also construct Goldstino at finite momentum: 



I^G(q)) = KHo) 



by boosting to finite momentum: 



(15) 



(16) 



Since the statistics of a single fermion created by has 
no physical consequence, |'0G(q)) is just like ferromag- 
netic spin wave states for SU{2) bosons[ll]; they have 
exactly the same dispersion and differ only in statistics. 
While IV'G(q)) is not an exact eigen state of Ha for finite 
q, it approaches one in the long-wave length limit q — > 0, 
and has quadratic dispersion: Eq — Eq oc |qpjjj]. 

It is appropriate at this point to discuss the relation 
between the sharp fermionic collective mode we call Gold- 
stino here, and the Goldstino in high-energy context. In 
high-energy context, Goldstino refers to the Goldstone 
fermion arising from spontaneous breaking of the global 
supersymmetry; it is a Weyl spinor with spin- 1/2. In 
our non-relativistic model, the gapless Goldstone fermion 
mode is also the result of the spontaneous breaking of the 
global supersymmetry. In this sense, they are quite sim- 
ilar. In fact, if we had a two-component boson system 
instead of a Bose-Fermi mixture, we would have a pseu- 
dospin ferromagnet that breaks SU(2) symmetry, with a 
branch of gapless spin-wave mode; the spin- wave mode is 
the Goldstone boson associated with breaking of SU(2) 
symmetry. In the Bose-Fermi mixture, the second com- 
ponent is fermionic, and the SU(2) symmetry becomes 
supersymmetry that is generated by the fermionic oper- 
ator Q. The difference here is that the Goldstino is a 
spinless fermion with quadratic (instead of linear) dis- 
persion, due to absence of Lorentz symmetry. 

We now turn to the more interesting case in which 
there is a finite density of fermions in the ground state 
\ipo)- In order to sustain this we must have a higher 
chemical potential for the fermions, thus A/i — — 
/is > 0. A/i can be viewed as a chemical potential for 
Goldstino; in its presence the fermionic Goldstino modes 
are "filled up" to some Fermi wavevector kp- 

In the presence of A/i > 0, Hq is no longer supersym- 
metric; it instead has a nonzero commutation relation 
with the supersymmetry generators: 

[Q,Hg\ - -[Q,^JiFNF + ^iBNB] = -A/iQ; 

[Q\HG\ = -[Q\^iFNF + ^lBNB\^^^iQ^■ (17) 



From (HT]) we can immediately conclude the following: 
(i) QIV'o) is an exact excited state with excitation energy 
A/x; and (ii) Q'^\'4i{i) = 0, because if it were not null, it 
would be a state with negative excitation energy —A/i. 
We thus find even though in this case supersymmetry 
is explicitly broken by A/i, we still have a sharp zero 
momentum fermionic collective mode generated by Q, 
which is now gapped. The situation is somewhat similar 
to what happens to a ferromagnet in an external mag- 
netic field: the field breaks rotation symmetry and opens 
a spin wave gap, but the spin wave remains a sharp col- 
lective mode. This is a "hole-like" Goldstino mode since 
it is created by Q instead of Q^, which creates a hole 
in the occupied Goldstino Fermi sea. The analogous (to 
Eq. (|15p ) finite momentum states are Rq\^Q), which are 
expected to have downward quadratic dispersion of the 
form _Eq — £'o « A/i — a|qp. 

Again due to (|17p . all eigenstates of Hg except for 
vacuum come in pairs whose energies differ by A/i and 
fermion numbers differ by one. This is because the canon- 
ical Hamiltonian remains supersymmetric. 

We now turn the discussion to possible experimental 
detection of the Goldstino modes. Normally one would 
expect that these modes can only be detected in processes 
in which a boson is turned to a fermion or vice versa, as 
that is what or Q does. There is no such process that 
can be easily engineered in cold atom systems that we 
are aware of, except for possible co-tunneling processes 
in which a fermion leaves the system and a boson enters 
the system simultaneosuly, or vice versa. In the follow- 
ing we show that in the presence of a Bose condensate, 
the Goldstino mode contributes a finite spectral weight 
to the spectral function of single fermion Green's func- 
tion; as a result it can be detected through processes that 
involve a single fermion. Physically this is possible be- 
cause in the presence of a Bose condensate, the boson 
number is not fixed in the ground state; as a result a sin- 
gle fermion (hole) can grab a boson from the condensate 
and propagate as the Goldstino mode. To demonstrate 
this we calculate the overlap between the fermionic single 
hole state /q^olV'o) with the normalized Goldstino mode 
(l/V^)g|V'o): 




q=0 



(18) 



where Ng is the number of bosons in the condensate 

and n^^o = ('/'o|/q=o/q=o|V'o) 1- When N% is macro- 
scopic (as is the case in the presence of a condensate), 
Q.olV'o) = ^/iV^j-^o), and the first term in dominates 
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the second. We thus find the zero momentum fermion 
Green's function has finite weight on the Goldstino, and 
the weight is approximately Ng/N, proportional to the 
condensate density. As a result the fermion spectral func- 
tion A{q = 0, Lu) has a sharp ^-function peak at w = A/x. 
This is highly unusual as in electronic or other fermionic 
systems, one normally expects sharp (coherent) spectral 
peak for fermions with momentum near kp and corre- 
sponding energy near zero in a Fermi liquid, correspond- 
ing to Landau quasiparticles which are well defined only 
near the Fermi surface. The sharp, coherent fermion peak 
at zero momentum and finite energy we find here is a re- 
markable consequence of the combination of supersym- 
metry and Bose condensation, which is unique to such su- 
persymmetric Bose-Fermi mixtures. Another remarkable 
property is that this spectral peak remains sharp at finite 
temperature T as long as T is below the Bose condensa- 
tion temperature Tc. This is because p7)) guarantees that 
the Goldstino mode is sharp at any T, while for T < Tc 
the condensate density is nonzero, so the fermion spectral 
function has a weight (proportional to condensate den- 
sity) on the Goldstino mode. We are not aware of any 
fermionic mode (like Landau quasiparticles) that remain 
sharp at finite T in other contexts. The fermion spectral 
function at small but nonzero momentum q will also have 
finite weight on the finite q Goldstino mode for T < Tc- 
This will result in a spectral peak at uj = A/i — a|qp, 
whose width goes to zero as |q| 0. 

In electronic condensed matter systems, single electron 
spectral function can be measured using electron tunnel- 
ing or photo-emission. In cold atom systems we do not 
have equivalent methods (yet) . On the other hand alter- 
native ways to measure single particle Green's function 
or spectral function are currently being developed, like 
stimulated Raman spectroscopy ^3j- Once these meth- 
ods become available, the Goldstino mode with the spe- 
cific properties discussed above can be probed through 
the fermion spectral function, as long as there is a Bose 
condensate (which is what bosons tend to form, except at 
certain commensurate lattice filling). We are thus highly 
hopeful that supersymmctric Bose-Fermi mixtures can be 
studied and the Goldstino mode can be detected experi- 
mentally, and our predictions can be tested. 

We close by stating that we have studied possible real- 
ization of supersymmetry in cold atom Bose-Fermi mix- 
tures, in a non-relativistic setting. Our study gives ex- 
amples of supersymmetry breaking, both spontaneous 
and explicit. In both cases the system supports a sharp, 
collective fermionic excitation which is a Goldstino-like 
mode; we have discussed ways to detect it experimen- 
tally as an unambiguous signature of supersymmetry. We 
hope such a study can provide hints to the mechanism 
of supersymmetry breaking in relativistic quantum field 
theories, in which the supersymmetry algebra is usually a 
graded Poincare algebra because of the linear dispersion 
of the relativistic fermion and that the boson is of the 



Klein-Gordon type. As a result the anti-commutation 
relation between the super-charge and its hermitian con- 
jugate is the Hamiltonian in that case, instead of the 
total particle number in the present work. 
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